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J: 

We prove the existence of axially symmetric solutions to the Vlasov- 
Poisson system in a rotating setting for sufficiently small angular velocity. 
The constructed steady states depend on Jacobi's integral and the proof 
relies on an implicit function theorem for operators. 
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CO ; 1 Introduction 

CO 

' In stellar dynamics, the evolution of a large ensemble of particles (e.g. stars) 

which interact only by their self-consistent, self-generated gravitational field, is 
CO ; described by the Vlasov-Poisson system 

OO 
O 



d t f+v-V x f-V x U-V v f = 0, (1.1) 
AJ7 = 4tt / 9, (1.2) 
p(t,x)= / f(t,x,v)dv. (1.3) 



Here / = f(t,x,v) > is the phase-space density, where £€R denotes time, and 
x,v£M. 3 denote position and velocity. U = U(t,x) is the gravitational poten- 
tial of the ensemble, and p = p(t,x) is its spatial density. We are looking for 
stationary solutions of Ql.l[) - (|1.3[1 . The ansatz 

fo(x,v) = <i>(E) = <P(^v 2 + U(x)) (1.4) 

is well known and autmoatically satisfies the Vlasov equation (|l.lj) , because the 
particle energy 

E{x,v):=\v 2 + U{x) 
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is a conserved quantity along characteristics. But we still have to construct 
the self-consistent potential. This is done by plugging (|1.4[) into the Poisson 
equation, more precisely, we have to solve 

A?7 = 471-/1$ ([/) = 4tt J H^v 2 + U(x))dv (1.5) 

The ansatz (|1.4j) only leads to spherically symmetric stationary solutions 
of (|1 . ![) - (|1 .3[) . where / is called spherically symmetric, iff f(Ax,Av) = 
f{x,v) Vi£0(3). Indeed, this is a special case of a more general result of 
Gidas, Ni and Nirenberg, cf. [2|. If one is interested in stationary solutions with 
less symmetry, more invariants can be added to (jl.4p . so that the right-hand 
side of (|1.5|) explicitly depends on x. 

One possibility is to consider a rotating system. If the ensemble is rotating 
around a given axis, say the X3-axis, we can change to the rotating frame and 
change coordinates as follows: 

(:=R t x, Tj:=R t v-nx (R t x), 

where 

(cos(wi) sin(wi) 0\ 
— sin(wi) cos(wi) , fi:= 
l) 

and the (rotational) velocity uj > is given. The Vlasov-Poisson system then 
takes the form 

d t f + r 1 -V c f-(V c U + nx(nx() + 2(nxr ] ))-V v f = 0, (1.6) 
A c U(t,C) = ^p(t,0, (1.7) 

p(t,0 = jf(t,(,v)dv (1-8) 
and the characteristic system of the Vlasov equation (|1.6|) reads 

C = v 

r) = -d c U(t, C)- 2fix7? -fix (fix C) ' 

which has the following expression as a conserved quantity, if U is time- 
independent: 

where Ej is also called Jacobi's integral. A natural ansatz for the construction 
of stationary solutions of (|1.6p - (|1.8[) is now 

/(Cry) = <p(Ej) = y{\ \ V \ 2 + U(() \^r 2 ) (1.9) 



2 



for a suitable function (^:R^M + , where r := r(x) = vCi + Cf • I n the original 
coordinates x,v one easily verfies that this ansatz leads to 

g(x,v) := f(C,v) = <P(\v 2 + U{R t x)-ujP), 

where we define P as the third component of the angular momentum, that is 
P :— x\V2 — X2V1, which is a conserved quantity of the characteristic system of 
the Vlasov equation (jl.ip . if U is axially symmetric with respect to the X3-axis. 
Obviously, the function f = f(C,v) then automatically satisfies (|1.6[) and one 
has to solve the Poisson equation, where we relabel C and r\ to x and v, 

AU = J ip{^v 2 + U(x)-^LU 2 r 2 )dv=:h{LJ,r(x),U(x)). (1.10) 

So if we construct an axially symmetric U solving (jl.lOp , the corresponding 
functions (g,U), with g defined as above also will be a stationary solution of 
(jl.ip - (jl.3p . Clearly, our ansatz for / satisfies (|1.6[) without any symmetry as- 
sumptions on U and this gives hope for the construction of stationary solutions 
with less symmetry, for example triaxial systems. 

Equation (jl.lOp has been studied, among others, by Vandervoort, cf. [5]. He 
observed numerically, that if tp is of the form 

l p(Ej) = (E -E J f + - 3/2 ) (1.11) 

then for 0.5 <(3< 0.808 there are triaxial solutions to (|1.10[) for sufficiently 
large ui. For small uj or j3> 0.808, all numerically constructed solutions are 
axially symmetric. Consequently, (|1.10|) seems to be of particular interest for 
the construction of ellipsoidal systems, but to our knowledge no self-consistent 
ellipsoidal systems to (|1.1[) - (11.3[) or (|1.6p - (|1.8p have been constructed analyti- 
cally yet. 

We will prove that there exist axially symmetric solutions to (jl.lOp for 
small u> under suitable assumptions on ip, where we treat the case /3 > 5/2 
in p. lip . For this purpose, we require, that for lu = 0, we have a nontrivial, 
spherically symmetric solution (fo,Ua) of (jl.lOp . Note, that in this case the 
righthand-side of (jl.lOp only depends on Uq. For oj^O, we want to apply an 
implicit function theorem to get solutions, which arise by deforming Uq, where 
certain symmetries are conserved. The central idea, which makes this approach 
work is to look for a solution U u as a deformation of Uq, i.e., = Uo(g(x)) 
for some diffeomorphism j on I 3 , and to formulate the problem in terms of 
finding zeros of a suitable operator T over the space of such deformations 
instead of the space of the potentials. Whereas the original problem (jl.lOp 
had to be solved in IR 3 , we will only need to know the deformation on a 
compact neighbourhood of the support of the original solution {foiPo,Uo), 
and this provides useful compactness properties. Furthermore, finite radius 
and finite mass of the constructed solutions then are just consequences of the 
corresponding properties of (fo,po,Uo)- 

Although the allowed perturbations for the potential Uq only have mirror 
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symmetry which would match a triaxial system, we have up to now no method 
to exclude axial symmetry with respect to the a^-axis for the perturbations 
constructed by the implicit function theorem. 

The approach described above has been used by Lichtenstein for proving the 
existence of slowly rotating Newtonian stars, as described by selfgravitating 
fluid balls, cf. [H [S]. A translation of Lichtenstein's approach into modern 
mathematical language is due to Heilig, cf. [3]. 

The investigations made there were applied to the Vlasov-Poisson system in 
[8], where stationary solutions to (|l.ip ~ (|1.3p of the form f{x 1 v) — p{E)\p{ujP) 
were constructed. There, the potential U a-priori was axially symmetric, so 
that the expression P = x\V2 — x 2 v\ is a conserved quantity with respect to 
the characteristic system. The procedure described there is the basis of our 
approach. 

This paper is organized as follows: In the next section we rewrite the problem 
in terms of finding zeros of the operator T, we then state the main result and 
prove it using an implicit function theorem. For this, we need certain properties 
of T which can be proved as in [5] , except some minor technical modifications 
and one lemma, where the symmetry of the allowed perturbations enters in. In 
Section [31 we generalize this important lemma dealing with properties of the 
operator c^T(0,0) to mirror symmetry. 



2 The main result 

The mappings, which leave our solutions invariant, are in the set 

S :={nio : (xi,sc2,ar3)i->(a:i,a;2,-a;3)) Tun : {xi,x 2 ,x 3 ) ^ (xi, -x 2 ,x 3 ), 
Ton : (xi,x 2 ,x 3 ) 1 * {-Xi,x 2 ,x 3 )}. 

Now let B R := {x G R 3 1 \x\ < R} and define 

C s (B R ):={fe C{B R )\ f(Ax) = f(x), AeS, xeB R }. (2.12) 

Then we have 

V/(0) = 0, if feC\B R )nC s (B R ). 
For tp:WL— > [0,oo[ we require 

{(pi) ipEC 1 (R) and there is i? ei with ip{Ej) = for Ej > E and <p{Ej) > 
for Ej<E Q . 

{(p2) ip is strictly decreasing in ] — oo,E [. 

{(p3) The ansatz fo{x,v) = p{Ej) with w = produces a nontrivial, spherically 
symmetric solution {fo,po,Uo) of (|1.1|) - (|1.3|) with po<=Co(R 3 ), supp po = 
B x and U eC 2 {R 3 ) with lim| a .|_ >00 t/'o(a;)=0. 



4 



Examples for a functions satisfying (ipl)-(ip3) are the so-called polytropes 

<p{Ej):=(Eq-Ej)X 
for k > 1 and suitable Eq < 0. Now we can state the main theorem. 

Theorem 2.1. Let r:= \Jx\ J rx\. There exists wo>0, such that for all coG 
]— loq,luq[ there exits a nontrivial solution (/ w ,U U ) of lll.6\) - fT78\) with 



(i) r(x,v)= 



<p(lv 2 + U u {x)-luj 2 r 2 ) for \x\<4 
else 



(ii) (f°,p ,U ) = (fo,Po,U ) and for \u\<w 0) (f a ,p w ,U") has the following 
symmetry properties: For all A£ S we have 

r(Ax,Av) = r(x,v), p"(Ax) = p(x), U^(Ax)^U^(x) 

and (f UJ ,U U ) is not spherically symmetric for lo =/= . 

(Hi) p^eC^R 3 ) and [/ u eC t 2 (R 3 ), where p w {x) = J f w (x,v)dv. 

(iv) The mappings ]—luo,ujo[3llh—> p u and ] — loo,ujo[3 u<— ► are continuous 
with respect to the norms \\ ■ ||i l00 or \\ ■ \\2,oo, respectively. 

Remark. If we add rotations about the X3-axis to the set S, the proof of 
Theorem 12.11 still holds - we can essentially follow the proof given here, and 
this shows that the constructed solutions in Theorem 12.11 have to be axially 
symmetric a-posteriori. This follows by the uniqueness of the mapping given 
by the implicit function theorem, cf. Theorem 14. II 



For the proof of Theorem 12.11 we need some lemmata. 

Lemma 2.2. The spherically symmetric solution {fo,poyUo) has the following 
properties. 

(a) The potential Uq is given by 

U (x) = - f f^rdy = -pr [ ' s 2 p (s)ds-4Tr [ sp (s)ds, xeR 3 . 



\x-y\ ' \x\ 



|.r| 



(b) Pq is decreasing with po(0) > 0, Uq (0) > and for every R>0 there exists 
C>0 ; such that U Q (r)>Cr, r€[0,R], and U (l) = E Q . 

(c) p' is Holder continuous and /7oGC 2 (K 3 ), where M 3 :=M 3 \{0}. 
Proof. The formula 

Uo{r) = A ^o s Y S)dS 
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easily follows from the Poisson equation with spherical symmetry and since we 
require lim| a; |_ >00 Uq(x) = 0, the representation for Uq holds by uniqueness. As 
to (b), for u = we have fo(x,v) = fo(E) = fo{\v 2 + Uq(x)) and this implies 

Po(x)= I f {x,v)dv = h (U (x)):=4:wV2 [ ° <p(E)^/E-U (x)dE, (2.13) 

JM 3 Ju (x) 

where the function h is continuously differentiable and with ((pi), (</?2) we have 
h'(s)<0 for ,s<_E . Consequently, p a is decreasing because U is increasing 
and since the steady state (fo,U ) is assumed to be nontrivial, we must have 
p (0) > 0. Thus actually U (r) > 0, r > 0, and since U o '(0) = (4tt/3)p (0) > this 
implies the estimate on Uq from below. The assertion that Uo(1) — Eq follows 
from (|2.13p and the assumption supppo = B\. The regularity of U follows from 
the formula for Uq above and the fact that pqEC^, which we deduce again from 
(|2.13| . Finally, the Holder continuity of p' will be part of the next Lemma. □ 

Lemma 2.3. Let E\ :=Uq(2) — Eg and define f by 



(p(±v 2 + U(x)-±L0 2 r 2 ) for U(x)<E Q + E 1 
else 



where tp satisfies {(pi), (ip2) and UeC 2 (R 3 ) with U(x)>Eq + Ei for |i|>4. 
Then the following holds: 

Pf( x ) : = / f{x,v)dv 
Jr 3 



with 



-h(cu,r(x),U(x)) 

{h{U{x)-\uj 2 r 2 ) for U{x)<E + E 1 
1 else 

^E-s(p(E)dE. 



(2.14) 



Furthermore, h€.C 1 (M.x [0,oo[xR) and for every bounded sef5cRx [0,oo[xR 
there are constants C>0 and /x€]0,l[ such that for (u>,r,u), (u)',r,u') G B we 
have 

\d r h(ui,r,u)\ < Cr, 

\h{uj,r,u)-h(Lu',r,u')\<C(\oj-Lu'\r + \u~u'\), 
\d u h(uj,r,u)-d u h(uj',r,u')\ < C(\u-u/\ + \u-uY)- 

In addition, for lo = 0, the function h(0,-,-) does not depend on r(x) and we can 
write ho := h(0,0,u). 



G 



Proof. Introducing polar coordinates, we have for U(x) < Eq + E 
p(x) = J ^(^v 2 + U(x)-^uj 2 r 2 ^dv 

= 4tt^ t 2 <p(^t 2 + U(x)-^u) 2 r 2 ^dt 
= 4ttV2 I (e 



E / i \l/2 

U{x) + -uj 2 r 2 \ <p(E)dE, 



and (gHH) follows. 

We have heC l CM) with 



f Eo 1 

for s < Eq and /i'(s) = for s > Eq and the first two estimates follow. Next, 

j pE -s y 

h"(s) = -4W2— / — —(p(E + s)dE 
dsj 2-jE 

= -47rV2/ — ^^'(S + s)^ 
Jo 2V-E 



yields local Lipschitz continuity of d u h with respect to w and u and the proof 
is complete. □ 

We want to find solutions of the equation 

AU = A7rh(L0,r(x),U) (2.15) 

and the main idea is to rewrite problem (|2.15p in terms of finding zeros of 
an operator T, which does not act directly on the space of potentials, but on 
deformations of the given spherically symmetric potential Uq- We define Banach 
spaces, which will serve as domain and range of T 

X:={/eC s (i34)|/(0)-0,/eC 1 (S4),3C>0:|V/(x)|<C7, a; ei?4, 

VxEdBi: lim ^\/f(tx)=:Wf(Qx) exists, uniformly in xEdBi}, 

where dB\ :— {x S M 3 | |x| = 1} and B± := £?4\{0}. We equip X with the norm 

11/11*:= sup |V/(aO|, feX 

and 

Y := {/ E C s (B 4 )\f(0) = Q,fe C\B A )^C > : |V/(s)| < C\x\,xeB 4 , 

VfCtx) V/(0x) 
Vx € dBi : lim =: exists, uniformly in x £ dB\\ 
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with norm 

ll/l|y: = s „ P ^fM, feY. 

To state more precisely, how to use functions in X to deform the potential Uo, 
we need the next lemma. 

Lemma 2.4. For(eX let 

5c :S 4 ^K 3 , g c (x):=x + ((x)-^- xeB 4 , Sc(°) = ° 

\x\ 

Then there exists r>0, such that for all (£^l, where 

n:={(eX\\\(\\ x <r} 

we have: 

(a) gQ-.Bi^Bi^-.— g^Bi) is a homeomorphism, <7f : -B4 — » -B^ is a C 1 - 
diffeomorphism, with 

\Dg c (x)-id\<^, xeB 4 

and for every x G dBi the mapping 

gC-OMBy^gciy) e 0, \g<;{4x)\x 

is one-to-one, onto and preserves the natural ordering of points in 0,4a;, 
where we defined xi,x 2 := {x\ + X(x 2 — x\) | A <G [0, 1] } for ii,i 2 €R 3 . 

(b) ||a;|<|5 C (a;)|<||a;|,a;eS4 } and g c (B 2 ) C B 3 , B 3 Cg c (B 4 ) C B 5 

(c) g i (Ax) = Ag i (x),xeB i and g^ 1 (Ax)^Ag^ 1 (x),xeB iX ,AeS 

(d) \Dg' ( 7 1 (x) — id\ < |, x € B^q and there exists a constant C ' > , such that for 
all (X' eft: 

± \g< (x) - g c (x) | + 1 Dg c (x) - Dg c (x) \ < C\\ C - C' II x , xeB 4 , 
\x\ 

and 

\g^(x)-g^(x)\<C\\(-C\\x\x\, xeB 3 
Proof In B4, we have for i,j= 1,2,3: 

d Xi gM = 5 ij+ d Xi «x)j± + ^ - ^f) (2.16) 

and therefore 

\Dg c (x)-id\<3\\C\\x. 
8 



With the inverse function theorem the first two assertions in (a) follow. For 
xedBi, 

g^tx) = tx + C,(tx)x = x(t + (,(tx)) 

and 

— (i + C(te)) = l + VC(te)-x>0 for || C||x small 
and the proof of (a) is complete. 

We have \(,(x)\ < \\(\\x \x\ for and this implies (b) for r>0 sufficiently 

small. Assertion (c) is easily verified, too. If we choose r even smaller we also 
have the first claim of (d), because 

Dg- 1 (x) = (Dg c r 1 (g^(x)). 

The estimate for g^ — g^ follows from the definition of g^ and the estimate for 
Dg c -Dg c follows from 

For xQBs, we have with (b): x £ g^(B 4 ) Pig^ (B4). Consequently, there exists 
y€B 4 mit x = g^{y). Now we have 

\g^(x)-g^(x)\ = \g^(g cl (y))-y\ 

= \9c 1 (9c(y))-9c 1 (ady))\ 

<2\g c (y)-g c (y)\<2U-C\\x\y\ 

<4||C-C'llxN, 



where we used the mean value theorem, the estimate for Dg^ and g^ (y ) , g^ (y ) C 
g c (B 4 ). ' □ 

We want to find solutions of (|2.15[) with the following structure 

U(x) = U ( (x) := Uo(g^(x)), xeB 4X , 

with a suitable (e£1. Obviously, we need U on the whole space M 3 , but this 
is only a technical problem. We use the fundamental solution of the Poisson 
equation to integrate (|2.15p and we then have to solve 

U (x)+ - <• dy = 0, xeB 4 . 2.17 

Jb 4i( \9d x )-y\ 

This equation essentially contains the operator we are looking for, but we have 
to modify things a little and also we want to get rid of the dependence on ( in 
the integration domain. 



Proof of Theorem \2.1\ For £ £ fl and we define 
T{w,Q{x):=U (x)+ — — — ^ dy 

Jb 3 \9d x )-y\ 
urn ( ~ h ^Ay),M9c\y))) , cR , 91jn 

-U (0)- / n — ^ dy, x£B 4 . (2.18) 
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Suppose we already know that this defines a continuous operator 

T-]-u>,u>[xil->Y 

for some Q > and T is continuously Frechet-differentiable with respect to £, 
where 

d c T(0,0):X^Y 

is an isomorphism - the first two assertion follow from [8] , Section 2 and the last 
assertion will be verified here in Section [3J It is also there that the symmetry 
of the perturbations plays a crucial role. 

The definition of Y requires T(w,£)(0) = and therefore we substracted the 
constant in (|2.18[) . With assumption (<^3), we know T(0,0)=0, because go = 
id and supppo = supp/io ° Uq = B\ C B 3 . The implicit function theorem, cf. pQ, 
Theorem 15.1, also stated in the Appendix as Theorem 14.11 cf. Section 2J now 
guarantees the existence of u>\ s]0,cD[ and the existence of a continuous mapping 

such that 

and C° = 0. We also will require that uj 2 r 2 <E\ in B4, where E\ is defined in 
Lemma |2"31 and therefore define 

w :=mm w b — - — \. (2.19) 

Now let £ = where we choose a fixed oje] — ujq,uiq[ and define 

p c (x):=h(u,r{x),Uo(gc l (x))), xeB 3 . (2.20) 

Then we have p$ € Cs(B 3 ) nC 1 ^). By Lemma [231 Pc >0 at most, if 
C/o(^ 1 (a;)) <i?o + ^i, which is equivalent to (x)\ <2 by Lemma [2~2l Con- 
sequently, 

suppp c =.g c (S 2 )c J B 3 . 
We extend by to all of M 3 and we achieve 

p c eC c (R 3 ), su PP p c cS 3 . 

We want equation (|2.20[) to hold everywhere, but we have not defined globally. 
We can rewrite T(oj,() = as 

U (x) = - [ . dy + C, xeB i; 

Jb 3 \9d x )-y\ 

or 

U (g^(x))=- [ pQLdy + C, xeB 4 . c , 
; Jb 3 \x-y\ 
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where 



Now define 



C:=U (0) + [ P -M dy . 

Jb 3 \v\ 



\v-y\ 

Then we have U c G C^R 3 ) with 

U ( (x) = Uo{g^(x)), xeB 3 cB 4ti (2.21) 

and thus p c e C^R 3 ) and U c E C b 2 (R 3 ) with AU C = 4tt Pc in R 3 . 
Furthermore, 

AU c = 4nh{uj,r(x),U c (x)), xeB 3 cB 4X . (2.22) 
The last equation holds even in R 3 . We have to show 

p{(x)=h(w,r(x),U{(x)), xeM 3 , 
that is, U c (x)> Eq + Ex for x <EM. 3 \g c (B 2 ). We know 
AU c {x) = } xeR 3 \gc(B 2 ), 

limu|_ >00 Z7£(a:) = C and 

C/ c (x)=£ + £i, xedg c (B 2 ), 
U ( {x)>E + E 1 , x€B 3 \g c (B 2 ). 

Here we used (|2.2ip and the monotonicity of E/o(M) with Uq(2) = Eq + E\. If 
C <Eq + E\, we have a contradiction to the maximum principle. Therefore, 
C > E + Ei and again by the maximum principle: U^> E + Ei on R 3 \(^(_B 2 ) 
and consequently, (|2.22j) holds in R 3 . 



Now define p u :=P£, U w := Uq and 
/"(ar,*): 



^(±u 2 + U^(x) - ±uj 2 r 2 ), for [f w (x) < £ + #i 



I else 

= U(±v 2 + U»(x)-±oj 2 r 2 ), for|x|<4 

[0 else. 1 ' 1 

Now defined by (|2.23p solves the Vlasov equation (|1.6p because it is constant 
along characteristics. More precisely, we have Uq(x) — ^ui 2 r 2 > E in a neigh- 
bourhood of 9i?4, if we choose ojo sufficiently small as in (I2.19[) . If we then 
fix (x,v) with Ej(x,v) < Eq and consider a characteristic (X,V) going through 
(x,v) we conclude that if XGB4, we have X £ B4 for all time. On the other 
hand, if x ^ B4, we have X B4 for all time. 
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Altogether, assertions (i)-(iii) of the theorem follow, except the non-spherical 
symmetry in the case uj^O. Choose i£R 3 with p u (x) > 0, x\ := a^=0,X2 — x^ — 
0. Then there exists some n GM 3 , such that 

^r, 2 + U"(x)~^ 2 a 2 <E . 

Now if (f u ,U u ) were spherically symmetric, there would exist a rotation A 
around the X2-axis such that (Ax)i = (Ax)2 = 0andf UJ (Ax,Av) — f u (x,v). But 
the monotonicity of tp implies 

t(x,v) = V ( l -v 2 + U«{x) - \u?a 2 ) = <p{Ej[x,v)) 

+ <p(Ej(Ax,Av)) = if{\ 2 + U»(x)) = r(Ax,Av), 

which contradicts our assumption of spherical symmetry. With a similar 
argument, one can also show that the constructed solutions cannot be axially 
symmetric with respect to any axis in R 3 except for the X3-axis. Though our 
deformations only have mirror symmetry with respect to every coordinate 
plane, which would match a triaxial system, we would still have to prove that 
the constructed are not axially symmetric with repect to the X3-axis to 
construct triaxial solutions. 

The asserted continuity properties (iv) can be proved as follows: For 
xGBs we have 

\U u (x)-U"'{x)\<\ I U' Q 1 1 oo \g£ (x) - g£ (x) \ < C\\ C - \\x. 

By the implicit function theorem, £ w continuously depends on w with respect 

to the || • ||x-norm and we have ff(x) =h{bJ,r{x),U 0J {x)). 

Lemma T2.3I implies that p^ is continuous in u> with respect to || • ||oo and 

U"(x) = - [ pHLdy + U (0)+ [ ^dy, xeR 3 

Jb 3 \x-y\ Jb 3 \y\ 

implies the continuity of U w in uj with respect to || -||i,oo- Differentiating the 
above expression for yields the continuity of with respect to || • ||i ]0 o and 
therefore also the continuity of U u in the norm || • ||2.oo- d 

3 c^T(0,0) is an isomorphism 

In this section, we want to establish some of the assumptions needed for the 
implicit function theorem. We will prove the following result: 

Proposition 3.1. The mapping <9^T(0,0) :X—tY is a linear isomorphism. 
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Let u>2 '■= \f\E\\/k, where E\ is defined in Lemma |2~31 and let us recall from 
[8], Proposition 3.1 that the Frechet-derivative of T :] — cj2,^2[xf2 — > Y is given 

by 

[£> f T(u,,OA](z) = 

(r-A — | - n) d u h(u,r(y),udy))vu ( (y) ■ ^MrM9t l (y))dy 
\\9d x )-y\ \y\J \g c (y)\ 



B 3 



h{w,r{y),U Q {y))dy—A{x), x€B 4 , (3.24) 

\9d x )-y\ m 



wherect>e]-ct> 2 ,cJ 2 [,CeO,AeX, and U Q {y) := U {g c 1 (y)), yeB 3 

We abbreviate Lo A :=c\T '(0,0) A for Ael. We observe that go = id and there- 
fore the function in (|3.24|) coincides with the potential Uq of the spherically 
symmetric steady state we started with, if £ = 0. We have 



p'd\x\) = d u h(0,r(x),U o (\x\))U^\x\) 

= d u h(0,r(x),U (\x\))\7U (x)~, xe 

\x\ 



J>3 



This implies 

(L A)(x) = -[ f^—-J-) p ' Qy\)A(y)dy-[ Jp> po (\y\)dy~A(y) 

JB 3 \\x-y\ \y\J JB 3 \ x -y\ fI 

= -U' {\x\)K{x)- f C_J__ 1 l-V (| y |)A( ! /)di/ J xeB 4 ,AeX. 
JB 3 \\x~y\ \y\J 

Now let 

(M)(I) ^WI 3 (^^) P ° (M)AM * xeB 4 ,AeC s (B 4 ). 
Then we can write 

{L A)(x) = -Ud\x\)[(id-K)A](x), xeB^AeX. (3.25) 
In order to prove Proposition l3.il we need 

Lemma 3.2. The linear operator K :C's(B 4 )—>Cs(B 4 ) is compact, where 
Cs{B 4 ) is equipped with the supremum norm || • ||oo- 

Proof. For AeC s (B 4 ) let 

V A (x):=- [ —!—pd\y\)A(y)dy, xeR 3 . 

Jb 3 \x-y\ 
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Then V A € C 1 ^ 3 ), VVa(O) = 0, and 

1 



(KA)(x)- 



(V A {x)-V A (0)), x€B 4 . 



Using Lemma \2.2{ c), we obtain the estimate 



\(KA)(x)\ < 



C\x\ 



IVyAllooN^CHAHoo, X £ -B4, 



where the constant C depends on p and Uq, but not on A or x. Thus K maps 
bounded sets into bounded sets. We next show that KA is Holder continuous 
with exponent 1/2, uniformly on bounded sets in CsiB^). Let M > and assume 
|| Alloc < M. In the following, constants denoted by C depend on po,Uo and M, 
but not on A. Obviously, p' AGL°°(R 3 ) and we deduce from Lemma l4~2l the 
existence of C > with 

\VV A (x)~VV A {x')\<C\\p , a A\\ 00 \x-x'\ 1 / 2 , x,x'eB 4 

Since VVa(O) = 0, the latter implies 

|VVa(^)| ^Cl^l 1 / 2 , xeB 4 . 

Now let x,x' e B A and |x| < |x'|. Then 

1 1 



\(KA)(x)-(KA)(x')\< 



05(1*1) u' (\x>\) 



\V A (x)-V A (0)\ 



1 



+ W |yAW -" AMI=:/l + l2 
and we obtain for some z£B± with \z \ < \x'\ the estimates 

h < m ^-™^\vvmto < C\x-x^M±^l ]z ^ 

<C\x-x'\^ 2 , 



and 



/ 2 < £n\VV A (z)\ Ix-ar'l < < C |x-.t'| 1/2 , 



so that 
and 



\(K A )(x)K A )(x')\<C\x-xf/ 2 , x,x'eB 4 
\(K A )(x)\ <C\VV A (z)\ <CM 1/2 , xeB A . 



We have shown that K maps bounded sets of Cs{Bn) into bounded and equicon- 
tinuous subsets of Cs{Bi). Thus K is compact by the Arzela-Ascoli theorem 
and the proof is complete. □ 
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Lemma 3.3. id — K : C's(Bi) — > Cs(Bi) is one-to-one and onto. 

Proof. Since K is compact, it suffices to show that id— K is one-to-one. Let 
AsCs(i?4) with A — A'A = 0. Now A = can be shown by expanding A into 
spherical harmonics. For that purpose, let 

{Snj, neN, j = l,...,2n + l} 

be the orthonormal set of spherical harmonics introduced in the Appendix, 
cf. Section 2J where for n £ N, the functions S n j : dB\ — » M, j = 1,. ..,2n+ 1 are 
homogeneous polynomials of degree n. We define 

A ni (r):= / 5„ li (0A(r0dwc = 4 / S ni j{x/r)h.{x)dLJ a (3.26) 



and we use the expansion of the integral kernel — y| into spherical har- 
monics, cf. Lemma 14.31 and Lemma 14.41 For x,j/eIR 3 , a; = r£ and y = sr] with 
^,?7e9i?i, r,s€R + , r^s, we have 

1 oo 2n+l . / • / \ \ n 

F-2/I ^ ^ 2n+l Vmax(r,.s)y 

ATA — A = then implies 

4n 1 f 3 2 , min (r,s)™ f 

" ^o( s ) ; — ^TT / S nd [Ti)A{sr])du n ds 



2n+\U' (r) 7 ru v ' max(r, s)«+i 

4tt 1 f 3 2 0" 
+ 2n+l[/ (r),/ S ^ S V+i 
4tt 1 f 3 2 s ( min(r,s) Tl 0" 

where we used that the functions S n j are orthonormal with repsect to 
(•,-)l2(9Bi)- We find that 



s Pa( s )-^TT / 5„, j (r/)A(s77)do; 7) ds 
1 2 ,,J min(r,s)" 0" \ . , , , 

S p (s) \ . . , i —r ]Anj(s)ds, 

Fln J \m&x(r,s) n+1 s n+1 J JW 



4tt [ r 

A i(r) = — 77777/ Po(s)s(s-r)A i(s)ds 
ru ov ) Jo 

and we obviously have lim r _»o Aoi (r) = 0. Let i? > be maximal such that Aoi (r) 
vanishes on [0,R]. Then for re [R,3], 

4tt /"" 
|A i(r)|< J lPolloo sup |A i(s)| / s(r- s)ds<C(r- R) sup |A 01 (s)|. 

rU o( r ) 0<s<r JR 0<s<r 

Thus for small e > 0, we have A i (r) = on the interval [R, R + e] and we conclude 
that Aoi vanishes on the whole interval [0,3]. Now up to linear combinations, 
the spherical harmonics for n=l are given by x±,X2,X3, and As Cs implies 

/ ^A(rOdu e = - [ aA«)^ = 0, 

JOB! JdBx 
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where we made the transformation £ h- >• (— ^1,^2,^3)- Analoguously, 



and we have An = A12 = A13 = 0. Let n > 2. Then 

47r 1 / C s n f 3 r" \ 

An]{r) = ^2^Tiu^){j Q s >o(s)^+rA ni ( s )d s +| s 2 p o(s)^+rA n ,( s )d s J, 

and 

471-1 / l /" r s n_1 A 3 r "-i \ 



4tt 1 / 1 



^ 1 / 1 3, ,3 / , 



AnjUoo -»r (-po)W + -2 / s Po(s)rfs + rp (r) 



= 2^TT I|A "''' 1100 ' 

where we integrated by parts in the third line and used the fact that U'o{r) = 
73- fo s 2 po(s)ds in the last line, also recall from (|2.13|) that —p' (r)>0. 
Now 2n + 1 > 3 for n > 2 implies that A n j = for n > 2 as well and the com- 
pleteness of {S n .j} induces A = 0. We conclude that id— K is one-to-one as 
claimed. □ 

It is now clear that Lq : X — » Y is one-to-one as well - this follows from Eq. 
(|3.25p and the fact that t/o(r) >0 for r >0. So once we have proved the next 
lemma, the proof of Proposition 13.11 will be complete. 

Lemma 3.4. Lq-.X^Y is onto. 

Proof. Let g£Y and define q:=g/U^. We will show q£X. We have q£ 
C 1 (B 4 )nC s (B 4 ) and 



!V7 1^ |Vg(g)l . I / x. 



^■^f \Vg(x) \g(x)\ \ 



By definition of Y and since [To £ C 2 ([0,oo[) with [/q(0)>0 we have that for 
every x £ dBi , 

V«?(0x) 1 (/(Oa;)^ 1 



[^'(0) 2 uv ^'(0) 2 ' 
as i^0+, uniformly in x£ dB\. 
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Since X <zCs{B/C), there exists by Lemma T3.3I an element AG Cs{Bi) such 
that 

A-KA = -q=-^-. 

TP 

This implies that L A = g and thus that Lq is onto, provided Ae J. To see 
the latter we observe that A = K A — q is Holder continuous since KA is Holder 
continuous. If we now define V A as above in the proof of Lemma 13.21 we also 
conclude that V A G C 2 (K 3 ) and thus KA e C 1 ^). Denoting by H Va the Hessian 
of V A we obtain for each x € B4 a point zgO,,t such that 



|V(XA)(z)|< 



\Va(x)-V a (0)\+ ^ x ^ \VV A (x)\ 



U(,(\x\) 

<-^-\(H Va (z)x,x)\ + ^-\WV a (x)\<C\\D 2 V a \ 



Finally, for a; € <9i?i , we have 

V(*T A)(tx) = -^(^aM - ^a(O)) + — VV A (te) 



' (H Va (0)x,x)x+-^-D 2 V a (0)x, 



2Ug(0y VAW ' ' ^'(0) 

as t^0+, uniformly in xEdBx. We have shown that KAgX and this implies 
A = KA + qEX and the proof is complete. □ 



4 Appendix 

In this section, we firstly state the implicit function theorem which is used for the 
proof of Theorem 1 2. II Then we give a regularity result for the Poisson equation 
and finally introduce spherical harmonics and state two important lemmas: an 
addition theorem and the expansion of the integral kernel l/\x — y\ in spherical 
harmonics. 

Theorem 4.1. Let X,Y,Z be Banach spaces, UcX and V CY neighbour- 
hoods of XoG X and yo^Y respectively, F :U x V — > Z continuous and conti- 
nously Frechet-differentiable with respect to the second variable. Suppose also 
thatF{x Q ,y ) = and F£ l (x ,y )^Z,Y). 

Then there exist balls B r (xo) <zU , Bs(yo) <ZV and exactly one continuous map 
G : B r (xo) — * B$(yo) such that Gxq — yo and F(x,Gx) = on B r (xo). 

Proof. [J, Theorem 15.1. □ 
Lemma 4.2. Let n<p<oo and let p(x)£L p (W n ) with compact support. Define 

V P (x):=- I , p(y)dy 
Jm™ \x-y\ 
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Then for every < a < 1 — n/p we have V p £ C 1,a (M. n ) and 

|a^ p (a:)-9 i y p ( a; ')|<C(n,a,p)| 2 :'~ a; r||/|| ?) /:"(supp{p}) i ^-? 
Proof. [6 , Theorem 10.2. □ 
Some facts about spherical harmonics 

In the following, we use the notation of [7j and we will always consider 
the case, where the space dimension q is equal to 3. For rteN, consider a 
homogeneous polynomial H n of degree n, which satisfies 

AH n (x) = 0. 

Then for ^dB x := {x G R 3 1 \x\ = 1}, 

is called a spherical harmonic of order n. For each n, there exist 2n+l lin- 
early independent spherical harmonics, which we call S n> j, j = l,...2n+l, cf. 
[7], Lemma 4. We denote by {S n ,j, n = 0,...,oo, j = l,...,2n + l} the orthonor- 
mal set of all spherical harmonics, where we orthonormalize with respect to 
(■i ■) L 2 (dB!)- Then we have the following 

Lemma 4.3. For a fixed nSN and £, rye <9-Bi, we have 
2,1+1 9 +1 

3'=1 7F 

where P n (t) is the Legendre Polynomial of degree n. 

Lemma 4.4. Let x,y 6t 3 with x = y = rr], for suitable £,77 € dB\ andr,R& 
R. TTien we have for R> r 

1 00 



\x-y\ " ^W.^ 

and for R<r 



n=0 x 7 

where P n (t) is the Legendre Polynomial of degree n. 
Proofs can be found in [7], Theorem 2 and Lemma 19. 
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